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Abstract 
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I. INTRODUCTION 



The spin and pseudospin symmetry [1,2] observed originally almost 40 years ago as a 
mechanism to explain different aspects of the nuclear structure is one of the most interest- 
ing phenomena in the relativistic quantum mechanics. It plays a crucial role for a Dirac 
hamiltonian with realistic scalar S{r) and vector V(r) potentials, for nucleon spectrum in 
nuclei, for the existence of identical bands in superdeformed nuclei, etc [3]. The key feature 
of the pseudospin symmetry is based on the small energy difference between single-nucleon 
doublets with quantum numbers nr,i,j — i+1/2 and rir — + — £ + 3/2, where rir, i 
and j are the single nucleon radial, orbital and total angular quantum numbers, respectively. 
These quantum numbers are relabelled as pseudospin doublets; £ + 1 — £ is the "pseudo" 
orbital angular momentum, s = 1/2 is the "pseudo" spin and j = £±s is the total "pseudo" 
angular momentum for the two states in the doublet [4]. For example, "nrSi/2, (^r — 1)^^3/2" 
is valid for i = I, "nrP3/2, ijir — l)/5/2 is valid for £ = 2, etc. Another key feature is the 
single-particle Hamiltonian of the oscillator shell model. This means that the pseudospin 
concept in the nuclear theory is a division of the single-particle total angular momentum 
into pseudo rather than normal orbital and spin parts. The shell model implies that nucle- 
ons move in a relativistic mean field produced by the interactions between nucleons. The 
relativistic dynamics of nucleons moving in the relativistic mean field are described by using 
the Dirac equation and not the Schrodingcr equation. 

The pseudospin symmetry concept is investigated by the framework of the Dirac equation 
and occurs as a symmetry of the Dirac hamiltonian when an attractive scalar S{r) and a 
repulsive vector V{r) potentials near equal to each other in magnitude, but opposite in sign, 
i.e., S{r) ~ —V{r). On the other hand, the sum of the vector and scalar potentials in the 
Dirac equation is a constant, i.e., V{r) + S{r) = constant, for the solution of the pseudospin 
symmetry in nuclei. This condition has been found by Ginocchio [5] and applied to the case 
of the spherical harmonic oscillator [6] . Meng et al [7] showed that the pseudospin symmetry 
is exact under the condition of d{y{r) + S{f))/dr — 0. Lisboa et al. studied the generahzed 
harmonic oscillator for spin-1/2 particles by setting either E(r) = V{r) + S{r) = or 
A(r) = V{r) — S(r) = [8]. A necessary condition for occurrence of the pseudospin 
symmetry in nuclei is to consider the case S(r') = [3,5,9,10]. For more realistic nuclear 
systems, the quality of the pseudospin symmetry is increased in the framework of the single- 
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particle relativistic models and hence the competition between the pseudo-centrifugal barrier 
and the pseudospin-orbital potential is completed in the onset of pseudospin symmetry 
[11]. The Dirac equation with the pseudospin symmetry is solved numerically for nucleons 
which move independently in the relativistic mean field with external scalar and vector 
potentials [12,13]. In addition to the numerical solutions, some analytical solutions are 
also discussed for solving the Dirac equation for some realistic potentials [14-17] with the 
pseudospin symmetry. The analytical solutions show that under the condition of pseudospin 
symmetry, the exact solution of the Dirac equation gives the bound-state energy spectra and 
spinor wave functions [18-20]. 

The aim of this paper is to present an analytical bound state solutions of the Dirac 
equation for the Hulthen potential under the conditions of the exact pseudospin symme- 
try and exact spin symmetry. To obtain a general solution for all values of the pseu- 
dospin (spin) quantum numbers, the pseudospin (spin) symmetry and orbital dependency, 
pseudospin-orbit (spin-orbit) dependent coupling are included to the lower component of 
the Dirac equation as an integer quantum number. This component has the structure of the 
Schrodinger-like equations with the pseudo-centrifugal (spin-centrifugal) kinetic energy term 
and its solution is analyzed by using some algebraic methods and efi^ective approaches. One 
of these effective approaches is applied to the pseudo-centrifugal (spin-symmetry) kinetic en- 
ergy term in the case of ^ > (£ > 0) and also an effective potential suggested in the form of 
the square of the Hulthen potential is taken into account instead of the pseudo-centrifugal 
kinetic energy term. For small values of the radial coordinate r, this effective potential 
gives a centrifugal energy term in the first approximation. Therefore, the pseudo-centrifugal 
(spin-centrifugal) kinetic energy term is accepted as an effective term in this region. It is 
worthy to state that Jia et al [21,22] have proposed an improved new approximation scheme 
to deal with the centrifugal kinetic energy term in the solution of the Schrodinger-Hulthen 
problem. Using this approximation scheme, Jia et al [23,24] have obtained approximate 
analytical solutions for the Dirac-generalized Poschl- Teller and Klein-Gordon-Poschl- Teller 
problems including the centrifugal kinetic energy term. Recently, Ikhdair [25] has applied 
the approximation scheme to deal with the orbital centrifugal term in the Schrodinger- 
Manning-Rosen problem using the Nikiforov-Uvarov method. Further, the approximation 
has also been applied to the Schrodinger-Hulthen problem using the improved quantization 
rule [26]. 
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In the present work, the Dirac equation for the Hulthen potential is arranged under 
the condition of the exact pseudospin (spin) symmetry and it's solution is obtained sys- 
tematically by using the Nikiforov-Uvarov (NU) method [27-31]. As an application of the 
Dirac-Hulthen problem with the pseudospin (spin) symmetry, the relativistic eigenvalue 
spectrum for various degenerate states is presented for several pseudo-orbital (spin-orbital) 
and pseudospin (spin) quantum numbers. 

The structure of the paper is as follows. In Sec. 2, the basic ideas of the Nikiforov-Uvarov 
(NU) method are outlined in short. In Sec. 3, the Dirac equation is briefly introduced for the 
spin and pseudospin symmetry solutions. In Sec. 4, the Hulthen potential is substituted into 
the lower component of the Dirac equation and the pseudo-centrifugal (or spin-centrifugal) 
kinetic energy term is replaced by the square of the Hulthen potential to apply the Hulthen 
square approximation. The main results obtained in previous sections are connected by 
means of the main equation of the NU method. Lastly, the general procedures of the 
solution method are followed to obtain the energy eigenvalue equation and two-spinor wave 
functions. Results and conclusions are performed in Sec. 5. 

II. BASIC IDEAS OF THE NIKIFOROV-UVAROV (NU) METHOD 

It is especially well known that the solutions of the Schrodinger and Schrodinger-like 
equations including the centrifugal barrier and/or the spin-orbit coupling terms have not 
been obtained straightforwardly for the exponential-type potentials such as Morse, Hulthen, 
Woods-Saxon, etc [32]. Although the exact solution of the Schrodinger equation for the 
exponential-type potentials has been obtained for i = 0, any £-state solutions have been 
given approximately by using some analytical methods under a certain number of restrictions 
[33,34]. One of the calculational tools utilized in these studies is the NU method. This 
technique is based on solving the hypergeometric type second-order differential equations 
by means of the special orthogonal functions [35] . For a given potential, the Schrodinger or 
Schrodinger-like equations in spherical coordinates are reduced to a generalized equation of 
hypergeometric type with an appropriate coordinate transformation r — > s and then they 
are solved systematically to find the exact or particular solutions. The main equation which 
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is closely associated with the method is given in the following form [27] 

+ v,'(.) + ^^(.) = 0, (1) 

where a{s) and 5(s) are polynomials at most second-degree, r(s) is a first-degree polynomial 
and ip{s) is a function of the hypergeometric type. 

Let us now try to reduce Eq.(l) to a comprehensible form by taking ip{s) — (f){s)y{s) and 
choosing an appropriate function (f){s): 

At the first stage, Eq.(2) can be seen to be more complicated than the main equation, 
Eq.(l). To ensure the reasonable understanding, the coefficient of y'{s) is taken in the form 
t(s)/(7(s), where t(s) is a polynomial of degree at most one, i.e., 



and hence the most regular form is obtained as follows, 

0'(g) ^ 

where 



(4) 



As)-^l[r{s)-T{s)]. (5) 
The most useful demonstration of Eq. (5) is 

r(s) = r(s) + 27r(s). (6) 

The new parameter 7r(s) is a polynomial of degree at most one. In addition, the term 
4>" (s) / (j){s) which appears in the coefficient of y{s) in Eq.(2) is arranged as follows 

In this case, the coefficient of y{s) is transformed into a more suitable arrangement by taking 
the form in Eq.(4); 

£(£) £(£)?(£) a(£^^a(£^ 

where 

a{s) = a{s) + 7r^(s) + 7r(s)[f (s) - a\s)] + 7r'{s)a{s). (9) 
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Substituting the right-hand sides of Eq.(3) and Eq.(8) into Eq.(2), an equation of the same 
type as Eq.(l) is obtained as 

/M + ^.'W + ^.W^o. (10) 

As a consequence of the above algebraic transformations, the functional form of Eq.(l) 
is protected by following a systematic way. Therefore, the transformations allow us to 
replace the function of the hypergeometric type ip{s) by the substitution 0(s)y(s), where 
4>{s) satisfies Eq.(4) whit an arbitrary linear polynomial n{s). If the polynomial a{s) in 
Eq.(lO) is divisible by cr(s), i.e., 

a{s) = Xa{s), (11) 
where A is a constant, Eq.(lO) is reduced to an equation of hypergeometric type 

a(s)y" + r(s)y' + Xy^O, (12) 

and also its solution is given as a function of hypergeometric type [35]. To determine the 
polynomial 7r(s), Eq.(9) is compared with Eq.(ll) and then a quadratic equation for 7r(s) is 
obtained as follows, 

n^s) + n{s)[T{s) - a'{s)] + a{s) - ka{s), (13) 

where 

A; = A-7r'(s). (14) 
The solution of this quadratic equation for 7r(,s) yields the following equality 



, , a'{s)—T{s) , / /cr'(s) — r(s) V -/ \ r / \ 
n{s) = ^ \ — — ± W ^ \ — —] -a{s) + ka{s). (15) 



In order to obtain the possible solutions according to the plus and minus signs of Eq.(15), 
the parameter k within the square root sign must be known explicitly. To provide this 
requirement, the expression under the square root sign has to be the square of a polynomial, 
since 7r(s) is a polynomial of degree at most one. In this case, an equation of the quadratic 
form is available for the constant k. Setting the discriminant of this quadratic equal to 
zero, the constant k is determined clearly. After determining k, the polynomial 7i{s) is 
obtained from Eq.(15), and then r(s) and A are also obtained by using Eq.(5) and Eq.(14), 
respectively. 
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A common trend which is followed to generahze the solutions of Eq.(12) is to show that 
all the derivatives of functions of hypergeometric type are also of hypergeometric type. For 
this purpose, Eq.(12) is differentiated by using the representation Vi{s) = y'{s) 

a{s)v';{s) + ri{s)v[{s) + fiM^) = 0, (16) 

where Ti(s) = r(s) + o-'(s) and fi^ = X + t'(s). Ti(s) is a polynomial of degree at most 
one and Hi is independent of the variable s. It is clear that Eq.(16) is an equation of 
hypergeometric type again. By taking V2{s) — y"{s) as a new representation, the second 
derivation of Eq.(12) becomes 

a{s)v'^{s) + T2{s)v'^{s) + nMs) = 0, (17) 

where 

r2(s) = ri(s) + a'{s) = r(s) + 2a'(s), (18) 

1^2 = 1^1+ r[{s) = A + 2r'(s) + a"{s). (19) 

In a similar way, an equation of hypergeometric type for Vn{s) — y^'^\s) is constructed as a 
family of particular solutions of Eq.(12) corresponding to a given A; 

a{sK{s) + TnisKis) + ii^Vnis) = 0, (20) 

and here the general recurrence relations for r„(s) and /x„ arc found as follows, respectively, 

T„(s) = T(s)+n(7'(s), (21) 

/.„ = A + nr'(s) + ^^^^a"(5). (22) 
When ji^ — 0, Eq.(22) becomes as follows 

A = A„ = -nr'(s)-^^^^^a"(5), (n = 0, 1, 2, . . .) (23) 

and then Eq.(20) has a particular solution of the form 

which is the Rodrigues relation of degree n and p{s) is the weight function satisfying the 
differential equation 

[cr(r)p(r)]' = T{r)p{r). 

To obtain an eigenvalue solution through the NU method, the relationship between A and 
A„ must be set up by means of Eq.(14) and Eq.(23). 



III. DIRAC EQUATION 



In the relativistic description, the Dirac equation of a single-nucleon with the mass // 
moving in an attractive scalar potential 5'(r) and a repulsive vector potential Vlr) can be 
written as 

'a.cP + Pii^c' + S{r)) + V{r)\ ^.^.(r) = i?n.«^„.«(f), (24) 



where 




P^-ihV, 5= , /3= , (25) 




with a is the vector Pauli spin matrix and / is the identity matrix. P is the three momentum 
operators, a and (3 are the usual 4x4 Dirac matrices [36] , c is the velocity of light in vacuum 
and h is the Planck's constant divided by 27c. En^^ denotes the relativistic energy eigenvalues 
of the Dirac particle. For nuclei with spherical symmetry, S{f) and V{r) potentials in Eq.(24) 
represent only the radial coordinates, i.e., S{r) — S{r) and V{r) — V{r), where r is the 
magnitude of r. The spinor wave functions ipn^^i^ ^e written in the following form 

Fn.«(r)[y,(^,0)x±]f^ 
Gn..(r) [n(^,0)x±]; 




where Yi{9, (j)) {¥^{0, 4>)) and x± are the spin (pseudospin) spherical harmonic and spin wave 
function which are coupled to angular momentum j with projection m, respectively. Fn^^{r) 
and GnrKif) are the radial wave functions for the upper and lower components, respectively. 
The label k has two explanations; the aligned spin j = £+1/2 (si/2,P3/2; etc.) is valid for the 
case of K = — (j + 1/2) and then I = £+1, while the unaligned spin j = £—1/2 {pi/2, (^3/2) etc.) 
is valid for the case of «; = (j -|-l/2) and then i = i— 1. Thus, the quantum number k and the 
radial quantum number rir are sufficient to label the Dirac eigenstates. The Dirac equation 
given in Eq.(24) may be reduced to a set of two coupled ordinary differential equations (in 
units oi c — h — 1): 

^ + ^) F^^,{r) = (// + - Air))GnAr), (27) 



- 7) GnAr) = ()" - Er^^K + S(r))F„,«(r), 



(28) 
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where A(r) = V{r)—S{r) and E(r) = V{r) + S{r) are the difference and the sum potentials, 
respectively. By substituting 

^"'^'^^''^ ^ (A^-^n.. + S(r)) (d^ " 7) 
into Eq.(27), the following second order Schrodinger-like differential equation for Gn^K(?^) 
can be obtained as 

(29) 

where En,.^. 7^ when E(r) = (exact pseudospin symmetry). Further, a similar equation 
for F„^fj(r) can be obtained as follows 

(30) 

where En^^^ ^ —/i when A(r) = (exact spin symmetry). Under the condition of exact spin 
symmetry, {dA{r)/dr = 0, i.e., A(r) = C ^constant), Eq. (30) turns out to be 

{£-2 - ^^75^ - (/^ + - C)S(r) + - /.^ + C (/. - K..)) i^n.«(r) = 0, (31) 

where i{£+l) comes from k,{k, + 1) and i{i+l)/r'^ is the spin-centrifugal kinetic en- 
ergy term. On the other hand, under the condition of the exact pseudospin symmetry 
{dT,{r)/dr = 0, i.e., E(r) = C ^constant), Eq. (29) is reduced to the form 

- + (/^ - + C)A{r) + El, -f,'-C{i, + K.«) j G„.«(r) = 0, (32) 

where i{£ + 1) comes from k{k — 1) and £{£ + l)/r^ is the pseudo-centrifugal kinetic energy 
term. According to the original definition of the pseudo-orbital angular momentum, the 
cases £ — K — 1 and £ — —k are valid for k > and k < 0, respectively. Therefore, the 
degenerate states come into existence with the same £ but different k, generating pseudospin 
symmetry. Another important point which is necessary to be said on Eq.(32) is that the 
radial part of the spinor wave function ipn^,{r) must satisfy the boundary conditions that 
Gurni^)/^ becomes zero when r ^ 00, and Gn^^{r)/r is finite at r = 0. 
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IV. BOUND STATE SOLUTION BY MEANS OF THE NU METHOD 



A. Hulthen Square Approximation 

In this section, we shall involve the Hulthen potential to solve the Dirac equation given 
in Eq.(32), meaning that the potential A(r) is exponential in r and the pseudo-centrifugal 
kinetic energy term is quadratic in 1/r. The exponential potential in r is the famous Hulthen 
potential [37,38]; 

p-Sr 

where 5 is the screening parameter which is used for determining the range of the Hulthen 
potential. The parameter Aq represents 5Ze^, where Ze is the charge of the nucleon [39]. 
The intensity of the Hulthen potential is denoted by Aq under the condition of 5 > 0. This 
potential has been used in several branches of physics and its discrete and continuum states 
have been studied by a variety of techniques such as the algebraic perturbation calculations 
which are based upon the dynamical group structure S0(2,l) [40], the formahsm of super- 
symmetric quantum mechanics within the framework of the variational method [41], the 
supersymmetry and shape invariance property [42], the asymptotic iteration method [43,44] 
and the approach proposed by Bicdcnharn for the Dirac- Coulomb problem [45,46]. With this 
potential in place, Eq.(32) has to be solved numerically because the exponential behavior of 
A(r) is not compatible with the quadratic behavior of the pseudo-centrifugal kinetic energy 
term. However, Eq.(32) is analytically solvable only for the zero value of the pseudo-orbital 
angular momentum, i.e., i — {k — 1). In order to obtain more realistic results relating to 
the degenerate states, the Dirac equation should be solved for any ^-states. In one of the 
methods used for solving Eq.(32), Hulthen square approximation can be introduced as an 
effective approximation to the pseudo-centrifugal kinetic energy term in the case of £ > and 
small r. Following the original work of Filho et al [41] for this approximation, an effective 
potential term can be considered as follows 

(l-e-sr)2 - ([1 + 5^ + ...]_ 1)2 - ^2 ■ y ) 

The exponential numerator in Eq.(34) is expanded for small values of r and higher-order 
terms are ignored up to first-order term. Recently, the authors of [42,44,46] have been used 
a more efficient approximation than that of Eq.(34) instead of the pseudo-centrifugal kinetic 
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energy term /r'^ . This approximation has the advantage that it is only vahd for small 

values of 5 and I. Whereas the present approximation in Eq.(34) can also be used for small 
values of 5 and I .. 

When A(r) is taken as the Hulthen potential and the approximation of the centrifugal 
term as in Eq.(34), Eq.(32) yields 



{H- + C) _ + h^^^ - n -C[n + L^^^) 



GnAr) = 0, 
(35) 



where k — £ + 1 ior k > and k — —i ior k < and the wave function has to satisfy 
the boundary conditions, i.e., Gn^r = 0) = and Gn^i^{r — > oo) = 0. It is convenient to 
introduce the following variable and parameters: 

s = e-^"^, r e (0, oo) ^ s e [0, 1] (36) 

2 ^IrK - -Gn- GEn^i^ , , 

OJi — ^2 ' K'^°) 

A^^uj + i^l-iie+l), (39) 
Bi^2ujI + uI, (40) 
which allows us to rewrite Eq.(35) in the simple form 

where the finiteness of our solution requires that G^^s = 1) = for r ^ and Gn^s — 
0) = for r — > oo. The above equation can be solved by using a special solution method 
mentioned in Ref.[27] and following a short-cut procedure given in Section 2. First of all, 
before starting the procedure of the solution, Eq.(41) is compared with the hypergeometric 
type differential equation given in Eq.(l) and consequently Eq.(41) is solved analytically due 
to the fact that the solution is still subjected to a methodology using algebra and calculus. 
This part will be treated in the next subsection partially. 

B. Pseudospin Symmetry Solution 

By applying the basic ideas of Ref.[27] and imposing the theory of orthogonal functions 
which are known as a generalization of the Rodrigues formula [35], the comparison of the 
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differential equations in Eq.(41) and Eq.(l) gives us the following polynomials; 

t(s) = 1 — s, a{s) — s{l — s), a{s) — Ais'^ — Bis + col. 



(42) 



In the present case, if we want to substitute the polynomials given by Eq.(42) into Eq.(15), 
the following equality for the polynomial 7r(s) is obtained 



7r(s) = -| ± i ^(1 - 4^1 - 4A;)s2 + 4(Bi + k)s - icul 



(43) 



The expression under the square root of the above equation must be the square of a poly- 
nomial of first degree. This is possible only if its discriminant is zero and the constant 
parameter k can be determined from the condition that the expression under the square 
root has a double zero. Hence, k is obtained as = 2a;^ — Bi ± iui (2£+l). In that 
case, it can be written in the four possible forms of 7r(s); 



7r(.) = -i ± i (- 



2£ + 1 - 2iu;i 
2i+l + 2iuji 



- 2iuji^ , for k+ = 2a;f - Bi + iui (2I+ , 
s + 2iuji^ , for k-. = 2ujI - Bi - iui (2i + . 



(44) 

One of the four possible forms of 7r(s) must be chosen to obtain an eigenvalue equation. 
Therefore, its most suitable form can be estabhshed by 



7r(s) = iuji — ^iuji s, 



for A;_. The trick in this selection is to find the negative derivative of t(s) given in Eq.(6). 
Hence, t(s) and t'{s) are obtained as 



t(s) = 1 + 2iuji - {2iuji + 2i+ 3)s, t'(s) = -{2iuji + 2^+ 3) < 



(45) 



In this case, a new eigenvalue equation for the Dirac equation becomes 



A. 



+ 2n, 



+ 2nriuji, 



(46) 



where it is beneficial to invite the quantity 



-rir-T s 



Tir(n,.-1) ^// 



a"{s) in Eq.(23). An 



other eigenvalue equation is obtained from the equality A = A;_ + tt' in Eq.(14), 



+ 1) (l + 2icc;). 



(47) 
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In order to find an eigenvalue equation, the right-hand sides of Eq.(46) and Eq.(47) must 
be compared with each other. In this case the result obtained will depend on En^^ in the 
closed form: 

/ (l + 2n,) (i + l) +nl + iylY 
/ / X ■ (48) 

Substituting the terms of right-hand sides of Eqs.(37) and (38) into Eq.(48), the energy 
eigenvalue equation for En^^ can be immediately obtained; 

c+^^]E^^^+(^X-^^^-c^^^o, (49) 

where 




Y2 I 'l-rn \ Y 



[/ = (1 + 2nr) + 1 j + nl, (50) 

Y = 2(nr + £ + l), (51) 

T.t/+(£±^. (52) 

S 

The energy spectrum of the Dirac equation for A(r) = V{r) — S{r) — —Aqj^^^ is obtained 

by means of Eq.(49). In this case, the states with the same rir and i will be degenerate. The 
two energy solutions of the quadratic equation can be obtained as 



E^ ^ — 



6^2 AqT + CV) ± ^A5^{Ao{C + /i) - 5'T)(AoAf + 8^T)Y^ + 8\C + 2/^)2^4 

2(Af+W) ■ 

(53) 

For a given value of Ur and k (or £), the above equation provides two distinct positive and 
negative energy spectra related with or -B"^^, respectively. One of the distinct solutions 
is only valid to obtain the negative-energy bound states in the limit of the pseudospin 
symmetry. Before seeking the acceptable solution, it is useful to present some analogy 
about the energy spectra. 

Now, we are going to find the corresponding wave functions for the present potential 
model. Firstly, we calculate the weight function defined as [47-51] 

1 If t(s) . \ 2- 2i+i 



p{s) = exp (^y -^d. ) = (1 - sY'^' , (54) 



and the first part of the wave function in Eq. (4) : 



(t>{s) = exp ( / "^ds ] = s'^' (1 - sf^^ . (55) 
13 



Hence, the second part of the wave function which is the solution of Eq.(20) can be obtained 
by means of the so called Rodrigues representation 



ynM 



(2£+l) 



nr+2iLij\ 



2s), [0,1], 



1-.) 



(56) 



where the Jacobi polynomial Pn^^'^\x) is defined for Re(i/) > —1 and Re(/x) > —1 for 
the argument x G [— and c„^k is the normalization constant. By using GnrK{,s) = 
(p{s)ynr{s), in this way we may write the lower-spinor wave function in the following fashion 

Gnr^ir) = Cn^^ (exp(-ia;i5r)) (1 - exp(-5r))^^^ pi^"^''^^""^^-* (1 - 2exp(-5r)) 
= Cr,.^K^^—^—7--^ (exp(-icc;i5r)) (1 - exp(-5r))^"'^^ 



X ^— rir, rir + 2 ^iwi + £ + 1^ ; 1 + 2ia'i; exp(— (5r)^ , 



where 



The hypergeometric series 2F1 {^—rir, nr + 2 (^icoi +£+1^ ; 1 + 2ia;i; exp(— 5r)^ 



(57) 



(58) 



terminates 

for rir = and thus converges for all values of real parameters a;i > and £ > 0. In case if 
C = 0, then iuJi6 = ^/ {jJ. + En^^) (/U — En^^) with the following restriction -E„^k < /i required 
to obtain bound state (real) sohitions for both positive and negative solutions of En^f^ in 
Eq. (53). Now, before presenting the corresponding upper-component Fn^^{r), let us recall 
a recurrence relation of hypergeometric function 



ds 



2F1 (a;&;c;s) = 2F1 (a + 1;6 + l;c+ l;s) , 



(59) 



with which the corresponding upper component Fn^i^{r) can be given by solving Eq. (28) as 
follows 



(exp(— icui^r)) (1 — exp(— 5r)) 
(a* - £^n.« + C) 



+ 1 



^ 5 exp(— 5r) 



(1 — exp(— 5r)) ' r 



X 2F1 i^—nr^rir + 2 (^ui 1^ ; 1 + 2ia;i; exp(— 5r) 

nr + 2 (iuj^ + £ + 1) 5 (exp(-5r))''"^+^ (1 - exp(-5r))^+^' 



11. 



(l + 2ia;i)(//-K.« + C) 



14 



X (^l-nr,nr + 2 (^iui + i+^^ ■,2{l + iLJi);exp{-6r)^ , (60) 

where E^^^ ^ /i when C — 0, exact pseudospin symmetry and d^^^ is the normahzation 
factor. 



C. Spin Symmetry Solution 

This symmetry arises from the near equahty in magnitude of an attractive scalar, S{r), 
and repulsive vector, V{f), relativistic mean field, S{r) ~ V{r) in which the nucleon move 
[48-51]. Therefore, we simply take the sum potential equal to the isotonic potential model, 
i.e., 

= -^oy^. (61) 

along with the approximation given by Eq. (34) to deal with the spin-orbit centrifugal term 
i{i + l)/r'^. In the last equation, the choice of E(r) = 2V{r) V{r) allows us to reduce the 
resulting solutions of the Dirac equation into their non-relativistic limits under appropriate 
choice of parameter transformations [51]. Therefore, the spin-symmetry Dirac equation (31) 
becomes 

- "^[iH'-^T + + - - - - c (" - ^»-)] } ^-('■) = 0- 

(62) 

where k, = ifoi k, > and k = 1) for k < 0. It is convenient to introduce the following 

new variable and parameters: 

s = e-^\ r G (0, oo) -> s G [0, 1] (63) 

vl = + (64) 
5 

^2 = — ^2 ' (65) 

A^^ul-vl-^{^ + l), (66) 

^2 = 2ul - v\, (67) 

which allow us to rewrite Eq.(62) in a more simple form as 

(ff' \-s d A-ii^- - B-ii^ ^ oJ\\ 
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where the finiteness of our solutions require that -F„^«;(l) = and F„^k(0) — > 0. We apply the 
NU method following the same steps of solution in previous section to obtain the expressions: 

f (s) = 1 - s, a{s) = s(l - s), a{s) = A2S^ - B2S + ul- (69) 

To avoid repition, the functions required by the method for 7r(s), k and t(s) can be estab- 
hshed as 

7r(s) = iuj2 - {iuj2 + ^ + 1) s, (70) 
k = 2u;l-B2- iu2 (2£ + 1) , (71) 

and 

t(s) = 1 + 2iu2 - {2iuj2 + 2i + 3)s, t'{s) = -{2iu2 + 2£ + 3) < . (72) 

respectively, with prime denotes the derivative with respect to s. Also, the parameters A 
and Xn take the forms: 

^nl + 2nr (-^ + 1) + 2n^ic^2 and X ^ ul - {£ + 1) {1 + 2iu2). (73) 

Using the basic condition A = A„ followed by simple algebra, we obtain 

,2\ 2 



V 2(n, + £ + i) y ' ' ' ' y ) 



and then the energy eigenvalue equation is immediately obtained 

\ IT'S [ ri I 2S'So\ „ / S5^ 



^^yis) J^"-"l^+^j^"^'^+lf J -/x^+c/. = o, (75) 

where 

W ={1 + 2nr) (^ + 1) + nl, (76) 

Z = 2(n, + ^ + l), (77) 

S = W+^^^^. (78) 
The two energy solutions of the quadratic equation (75) can be obtained as 



^"'-'^ 2(E2 + Z25') 



(79) 
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For a given value of rir and k (or i), the above equation provides two distinct positive and 
negative energy spectra related with or -B"^^, respectively. One of the distinct solutions 
is only valid to obtain the positive-energy bound states in the limit of the spin symmetry. 

In our calculations for the spin symmetry wave functions, we firstly find the weight 
function: 

p{s) = (1 - sf^' , (80) 
and from which the second part of the wave function by means of Rodrigues formula as 

yuM = an^^s-'^-^ (1 - s)-^''^'^ ^ [.-+^^- (1 - sr"-''"-'] 



^p(2..„2m)(i_2,), ,e[o,i], (81) 

where the Jacobi polynomial Pn^'^\x) is defined for Re(z/) > — 1 and Rc(/i) > —1 for the 
argument x G [—1, +1] and a„^K is the normalization constant. Further, the first part of the 
wave function is being calculated as 

0(s) = s'^'' (1 - s)^+' . (82) 

By using -F„^k(s) = 0(s)y„^(s), in the spin symmetry case, we may write down the upper- 
spinor wave function in the following fashion 

FnAr) = an^n {expi-tuj^Sr)) (1 - exp(-5r))'+^ Pi'r"''+'\l - 2 exp(-5r)) 

= anrK- ^"^ ^"^ {exp{-iu25r)) (1 - exp(-(5r))^^^ 

rirl 

X 2F1 {-nr,nr + 2{iuj2+i+l);l + 2iuj2;exp{-5r)) , (83) 

where 



ioo2S = ^JC {E^^, + El^ > 0. (84) 

It is noted that the hypergeometric series 2F1 {—rir, nr + 2 {iu2 + + 1) ; 1 + 2iu2] exp(— 5r)) 
terminates for = and thus it converges for all values of real parameters c<;2 > and 



£ > 0. In case when C = 0, then iuoi5 = + -Ew^k) (/i — En^,^) with a restriction for real 
bound states that < /i for both positive and negative solutions of En^^^ in Eq. (79). 
Thus, the corresponding spin-symmetric lower-component Gn^K{f) can be found as follows 

(exp(— ja;2()>)) (1 - exp(-(5r))^'^^ [(£ + 1) (5exp(-5r) _ ^ , n 

^rirKy) — Oji^i^ 



(// + - C) 
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(1 — exp(— or)) r 



X {-rir, nr + 2 {iuj2 + ^ + 1) ; 1 + 2iuj2; exp{-Sr)) 



UrS K + 2 (£ + 1 + iu2)] (exp(-^r))''^^+' (1 - exp(-(5r)) 
(1 + 2iuji) ill + En^, - C) 



X 2-Fi ( 1 - n^, + 2 ( ia;2 + ^ + - J ; 2 (1 + ioj2) ; exp(-(5r) 



(85) 



where £^n^K; 7^ —A* when C = 0, exact spin symmetry and hn^^ is the normahzation constant. 

Let us finally remark that a careful inspection to our present spin-symmetric solution 
shows that it can can be easily recovered by knowing the relationship between the present 
set of parameters (a;|, z/|, ^2, -B2) and the previous set of parameters Ai^ Bi). This 

tells us that the positive energy solution for spin symmetry (negative energy solution for 
pseudospin symmetry) can be obtained directly from those of the negative energy solution 
for pseudospin symmetry (positive energy solution for spin symmetry) by performing the 
following replacements [48-51]: 

FnAr) ^ G^^,{r), V{r) ^ -V{r) (or Eq ^ -Aq), £(£ + 1)^ i{U 1) 

> Etn ^ -KrK, ^2 ^ ^1 and 4 ^ -ul (86) 

That is, with the above replacements, Eqs. (49) and (57) yield Eqs. (75) and (83) and vice 
versa is true. 

Let us now present the non-relativistic limit. This can be achieved when we set C = 0, 
Eo = (5 and using the mapping En^^ — A* — ^ En^^ and En^^ + — )> 2/^ in Eqs. (64), (65) and 
(74), then the resulting energy eigenvalues {in. h — c — e — 1 units) are 



[1 + 2nr) (£ + 1) 5 + nl5 - 2/1 



(87) 



2 (n, + £ + 1) 

Also, the wave functions in Eqs. (83) and (84) turns out to become 

RnA{r) = a^^,r-' (exp(- V-2/xE..,r)) (1 - exp(-5r))^+^ pi'^^^^^/''''+'\l-2 exp(-5r)) 



(2v/=2^/c5 + l)^^ -1 ( I r^T^ IS ^^^+1 

X 2F1 (-n„ nr + 2 (^^-2nEn,i/5 + £ + l) ; 1 + 2 V-2/xE„^^/5; exp(-5r)) , E^^t < 0. 

(88) 
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V. RESULTS AND CONCLUSIONS 



In the present study, the Dirac equation for the Hulthen potential is approximately solved 
under the condition of the exact spin and pseudospin symmetry within the framework of 
the relativistic mean field theory. By using the basic ideas of the NU method, the energy 
eigenvalue expression for the arbitrary pseudo-orbital angular momentum £ is obtained ap- 
proximately. The second-order differential equation given in Eq.(32) is solved by applying 
the Hulthen square approximation to deal with the pseudospin-orbit and spin-orbit cen- 
trifugal and kinetic energy terms l)/r^ and l)/r^. The energy spectrum for any I 
states is obtained analytically. Under the condition of the exact pseudospin and spin sym- 
metry limitations, the energy relations in the Dirac equation with equal scalar and vector 
Hulthen potentials are recovered to see degenerate states. 

The results obtained for this motivation show the orbital dependency of the Dirac equation 
for the Hulthen potential. Certainly, an analysis detailed by solving Dirac equation in rel- 
ativistic mean field theories needs to use a very large scale (~ 660 MeV) comparing to the 
nuclear physics scale (few MeV) ) in point of the intensity of potentials [6]. For this reason, 
the intensity of the potential, Aq, used in Eq.(33) is considered as 3.4 fm~^. The units 
ofh — c = e = l are used throughout the present work for the sake of simplicity. Hence, 
the energy eigenvalue expression given in Eq.(53) can be simply discussed by using a set of 
physical parameter values. In the below explanations, although the energy spectrums can 
be calculated in dimensionless or arbitrary units, the calculations are preferably made in 
fm~^ for the energy, mass, C and intensity of the potential. 

It is noted that the energy spectrum given by Eq.(53) indicates a family of the pseu- 
dospin symmetry Hulthen potential. Moreover, the analj^ical expression for Eq.(53) can be 
confronted with the results of [44] which is slightly in agreement with the result presented 
in Eq. (47) for the pseudospin symmetry solution. The results are only valid for small values 
of S and i (or k). This spectrum changes with the relevant quantum numbers as well as the 
screening parameter S. The variation of the energy spectrums {E^^,^ and E'^^^) according 
to the screening parameter S is shown in Fig. la and Fig. lb, with the choices of parameters 
C = —4.9 fm~^ and fi = 5.0 fm~^, which is in the range of nucleon mass value (~ 1 
GeV). Figure la indicates the positive-energy bound states, i.e., -E^^^? while Fig.lb shows 
the negative-energy bound states, i.e., E~^^. For a given value of and i, it is seen that an 
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increment on 5 leads to an increment on along the negative-energy direction whereas 
the same increment on 5 results with a reduction on E^^^^ along the positive-energy direction. 
The results presented in Fig. la show that the energy difference between the states is still 
small although the values of screening parameter 6 increases. Figure lb has two interesting 
results: The first one indicates that the negative-energy bond states appear with the large 
values of S and £. The reason of this aspect comes from the approximation mentioned in the 
previous sections. The second one belongs to the small values of S. For instance, E~^^ for 
lsi/2, 1(^7/2 and 1^^9/2 is vahd under the condition oi 6 > 0.09, 6 > 0.05 and S > 0.03, respec- 
tively. Therefore, E'^^^ still represents the negative-energy bound states for small values of 
6 when £ increases. These results can be also expanded on the other states of the Hulthen 
potential with the pseudospin symmetry. 

It is well-known that for the finite nuclei the constant C is adjusted to zero because 
each potential goes to zero at large distances. If the difference between the scalar S{r) and 
vector V{r) potentials equals to a given constant C, this is equivalent to adding the relevant 
constant to the relativistic energy and mass. The constant for the energy is unimportant 
because it does not affect the energy difference. Whereas the variation of C is equivalent 
to the variation of mass. This is more physically transparent for the pseudospin-orbit de- 
pendency of the Dirac equation. Moreover, in the case of the infinite nuclear matter, the 
constant C could be non-zero. The energy spectrum versus the mass /i is plotted by setting 
C = —4.9 fm~^, as shown in Fig.2. The variation of the energy spectrum for ^ = 1, ^ = 3, 
£ — 5 and ^ = 7 is presented by using 5 — 0.25. The radial quantum number is fixed to 
1 (rir = 1). In Fig.2, the dashed (red fines, see colour onfine) and sofid (blue lines, see 
colour online) lines represent -E^^^ E~^^, respectively. According to Fig.2, there are two 
different regions of energy spectrum versus the mass. For < /j, < 2.3 fm~^ the energy 
spectrum is in the negative region completely. With the £ increasing, E^^^ fan out along the 
positive part of the energy spectrum whereas E~^i^ overlaps in going from £ — 1 io £ — 1. 
The case of in the interval of < < 2.3 fm~^ represents the degenerate states in 
the different values of £ for a given value of n^. A similar trend is seen when 5 > 3.2 fm"^ 
in the several values of £ for E^^^ by crossing the zero axis toward the positive direction of 
the energy spectrum. However, the numerical results of -E^^^ over the axis are not relevant 
for the negative-energy bound states. Meanwhile, the results of E~^^ are also valid under 
the conditions ol 5 > 3.9, 3.5, 3.3, 3.2 fm'^ for £ — 1, 3, 5, 7, respectively, but only 
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relevant in the large values of i. Furthermore, the energy spectrum versus the constant C 
is plotted by taking // = 5.0 jwT^ and 5 = 0.25 as shown in Fig. 3. According to Fig.3, it is 
seen that the negative values of C show more strongly binding energies under the condition 
of C < —11.0 fm~^ for E^^^ (dashed lines) and El^^^ (blue line) in the whole values of £. 
Moreover, E^^^ (dashed line) still shows the negative-energy bound states on condition that 
— 8 < C < — 5 /m~^ up to the zero axis. In the case of E~^i^ (blue lines) the situation 
becomes different than that of -B^^^. The axis is crossed with changing i for -B"^^. The 
crossing points of axis from C — —8 fm~^ become large with the £ increasing. In Table 1, 
the degenerate states are presented for a Dirac particle within the Hulthen potential, with 
C = —4.9 fm~^, fi = 5.0 fm^^ and Aq = 3.4 fm^^. The several pseudo-orbital and radial 
quantum numbers are used in the numerical calculations to predict the orbital dependency 
of the Dirac equation under the condition of the exact pseudospin symmetry. As an example, 
the Dirac eigenstate lsi/2 with rir — 1 and k — —1 will have a partner which is denoted by 
the 0(^3/2 with — 1 = and k — 2. These states are called the pseudospin partner and 
degenerated with each other. 
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FIG. 1: The variation of the energy spectrum in units of fm ^ versus the screening parameter 6. 



FIG. 2: The variation of the energy spectrum versus the mass fi. All parameters are in units of 



FIG. 3: The variation of the energy spectrum versus the constant C. All parameters are in units 
of fm~^. 
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TABLE I: The negative-energy degenerate states in units of fm of the pseudospin-symmetry 
Hulthen potential for various values of n^, i and 6. For a special case, = 5 fm~^, Aq = 3.4 fm~^ 



c 


= — 


-4.9 fm 


-1 












i 


n, 


. S 


Degenerate States 




rir 


S 


Degenerate States 




1 


1 


0.025 


(lSl/2, 0^3/2) 


0.0963638 1 


2 


0.025 


(2si/2, lci3/2) 


0.0928939 






0.100 




0.0425738 




0.100 




-0.0103694 






0.175 




-0.0710009 




0.175 




-0.2174930 






0.250 




-0.2346580 




0.250 




-0.4920870 


2 


1 


0.025 


(lP3/2, 0/5/2) 


0.0912282 2 


2 


0.025 


(2p3/2, 1/5/2) 


0.0863238 






0.100 




-0.0363590 




0.100 




-0.1078600 






0.175 




-0.2930130 




0.175 




-0.4732160 






0.250 




-0.6351320 




0.250 




-0.9131390 


3 


1 


0.025 


(1^5/2, O97/2) 


0.0839128 3 


2 


0.025 


(2^5/2, I57/2) 
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